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Van der Waals Interaction in Uniaxial Anisotropic Media 

Pavel Kornilovitch 1 '^ 

Hewlett-Packard Company, Imaging and Printing Division, Corvallis, Oregon 97330 USA 
(Dated: 10 January 2012) 

Van der Waals interactions between flat surfaces in uniaxial anisotropic media are investigated in the nonre- 
tarded limit. The main focus is the effect of nonzero tilt between the optical axis and the surface normal on 
the strength of van der Waals attraction. General expressions for the van der Waals free energy are derived 
using the surface mode method and the transfer-matrix formalism. To facilitate numerical calculations a 
temperature-dependent three-band parameterization of the dielectric tensor of the liquid crystal 5CB is de- 
veloped. A solid slab immersed in a liquid crystal experiences a van der Waals torque that aligns the surface 
normal relative to the optical axis of the medium. The preferred orientation is different for different materials. 
Two solid slabs in close proximity experience a van der Waals attraction that is strongest for homeotropic 
alignment of the intervening liquid crystal for all the materials studied. The results have implications for the 
stability of colloids in liquid crystal hosts. 



I. INTRODUCTION 

One component of the interaction between colloidal 
particles suspended in a fluid is van der Waals (vdW), 
or dispersion, forces.— In the nonretarded limit (distances 
below approximately 1 fim) the dispersion forces are me- 
diated by longitudinal matter-like modes localized on the 
interfaces, while in the retarded limit (large distances) 
they are mediated by transverse light-like modes standing 
between the interfaces?- 2 Quantum-mechanical treatment 
of the nonretarded regime was pioneered by London 3 
and of the retarded regime by Casimir and Polder A Lif- 
shita^ introduced a macroscopic treatment of the problem 
where the electromagnetic properties of the bodies and 
the medium were characterized through dielectric per- 
mittivity functions. Dzyaloshinskii and Pitaevskii gen- 
eralized the Lifshitz theory to non-uniform bodies using 
the formalism of temperature Green's funstions. -~- - Later, 
van Kampen, Nijboer, and Schram developed a method 
of calculating the vdW interaction based on the knowl- 
edge of electromagnetic surface states, which simplified 
the treatment of complex geometries^ Those pioneering 
works started theoretical and experimental investigation 
of vdW forces, a research field that has remained active 
ever since ji2r— 

The present article concerns with vdW forces between 
nanosize particles in anisotropic media. It is motivated 
by experiments on colloids suspended in liquid crystals^ 
In the cited paper, clustering of plate-like clay particles 
in liquid crystal 5CB was observed by small-angle X-ray 
scattering. Apart from fundamental interest, suspensions 
of colloids in nonaqueous solutions are of increasing tech- 
nological importance, most notably as the work medium 
in the electrophoretic reflective displays. 

The dielectric anisotropy was introduced into the 
vdW problem, apparently for the first time, by Ki- 
hara and Honda. 15 Those authors considered electrically 



anisotropic solid particles separated by an isotropic liq- 
uid. Katsi£ generalized the temperature Green's func- 
tion formalism to anisotropic media. He considered 
two anisotropic particles separated by an isotropic liq- 
uid and calculated a vdW torque that rotates the two 
particles relative to each other. Kats also considered 
a cholesteric liquid crystal mediating vdW interaction 
between two isotropic bodies. Parsegian and Weiss 17 
studied two anisotropic bodies interacting across an 
anisotropic medium but with all three regions sharing 
one common optical axis. Smith and Ninham 1 " investi- 
gated the vdW forces applied by two anisotropic bodies 
on a film of twisted nematic squeezed between them. The 
torque exerted by the bodies was balanced by the elastic 
energy of the nematic, which determined the structure 
of the twist. Sarlah and Zumer 1 ^ considered two semi- 
infinite optically uniaxial bodies separated by a uniaxial 
medium with all three optical axes parallel. Those au- 
thors derived an analytical expression for the Hamakcr 
constant but did not consider vdW torques. 
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The main goal of the present paper is to investigate the 
effects of arbitrary orientation of the optical axis relative 
to particle surfaces. Due to the complexity of the general 
problem, treatment will be confined to the parallel-plate 
geometry in the nonretarded limit. In this case, the vdW 
energy and forces can be derived from the spectrum of 
longitudinal surface modes t 12 ' 13 The electrostatic model 
is formulated in Section |TT] and the general solution is 
constructed in Section IIVI For numerical evaluation of 
forces dielectric functions at imaginary frequencies are 
needed. Such functions for several materials are listed 
in Section [IIH In addition, a temperature-dependent di- 
electric model of 5CB is defined for both polarizations. 
The single-slab problem is solved in Section [V] and the 
two-slab problem in Section IVT1 Some implications of the 
obtained results and future directions are discussed in 
Section |VH] 
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FIG. 1. The model geometry. Two finite thickness slabs 2 
and 4 made of dielectrically isotropic materials are immersed 
in uniaxial liquid crystals 1, 3, and 5. The slabs are infinite 
in the x and y directions and parallel to each other. The 
optical axes of the liquid crystals are tilted away from the 
surface normal by respective angles #1,3, 5 and rotated around 
the z-axis by respective angles f/>i,3,5. 



II. MODEL AND METHOD 

The overall geometry is shown in Fig. [T] The system of 
interest consists of two parallel slabs 2 and 4 a distance 
L = L32 — L43 apart. All slabs' surfaces are parallel to 
the a:y-plane of the coordinate system, and the z-axis is 
perpendicular to the interfaces. The slabs are infinite in 
the x and y directions and have finite thicknesses in the z 
direction: c?2 and respectively. The materials of 2 and 
4 are assumed to be optically and dielectrically isotropic. 
They are characterized by the scalar dielectric functions 
£2(w) and £4(0;). [The formalism developed below could 
be easily generalized to anisotropic solid materials. How- 
ever, it would have obscured the main focus of the paper, 
which are the effects derived from the anisotropy of the 
medium. Such a generalization is left for future research.] 

The space between the slabs as well as outside the slabs 
is filled with dielectrically anisotropic media. The respec- 
tive regions are labeled r = 1, 3, and 5. The media are 
assumed to be optically uniaxial such as most liquid crys- 
tals. (The terms "uniaxial media" and "liquid crystals" 
will be used interchangeably throughout the paper.) The 
dielectric and optical properties are characterized by the 
parallel and perpendicular dielectric functions e T (u>,T) 
and e^-(oj,T). In uniaxial liquid crystals, both dielectric 
functions are strong functions of temperature T. The 
theoretical formalism allows the three materials to be 
different. However in all the examples considered below 
the material will be the same in all three regions. 

In each region r, the direction of the optical axis is de- 
fined by a polar angle r and azimuth angle ip r , as shown 
in Fig. [TJ Within Cartesian coordinate system, the di- 



electric tensor is defined as follows, (i) In the starting ori- 
entation the optical axis coincides with the z-axis. This 
implies e zz = e", e xx = e vv = e ± , and all off-diagonal 
terms are zero, (ii) The medium is rotated around axis y 
by an angle 0. (iii) The medium is rotated around axis z 
by an angle ip. As a result, the dielectric tensor in region 
r assumes the form 



~y* P vv P yz 
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All five regions are assumed to be spatially uniform, 
i.e. the dielectric functions are independent of x, y, and 
z. This implies no distortion of the director field and 
hence a zero bulk elastic energy. Physically, this cor- 
responds to either very weak or very strong anchoring. 
In the first case, the surface energy is zero, and the ori- 
entation is an arbitrary parameter. In the second case, 
the orientation is fixed and the vdW energy is a small 
correction to the dominant surface energy. In any case, 
the goal of the present calculation is to characterize the 
vdW energy for a given orientation which is considered 
a model parameter. 

In the surface mode method^iii - — the vdW energy 
is derived from the spectrum of electromagnetic surface 
modes. If W(q, lj) — is a spectral equation, then the 
vdW free energy is a sum over all quantum numbers q 
and bosonic imaginary frequencies £„ = {2-nkT / K)n 



F = kT Y, Y^MW{qMn)] 



(3) 



The prime at the sum sign indicates the n — term must 
be taken with weight 1/2. Differentiating with respect to 
the distance between the bodies yields the interaction 
force. 

The method is particularly simple in the non-retarded 
limit when surface modes are obtained from solutions of 
the anisotropic Laplace equation rather than of the full 
set of Maxwell's equations: 
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(4) 



The modes are fixed by appropriate boundary conditions 
at the phase boundaries. This is completely analogous to 
a problem in optics. Mathematically, the solution can be 
constructed by using the formalism of transfer matrices. 
This will be done in Section HVl 
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TABLE I. Temperature-dependent coefficients C of the three- 
band dielectric model, Eqs. (0-©, of 4-cyano-4-n-pentylbi- 
phenyl (5CB).— Other model parameters are: uio = 9.19 eV, 
wi = 5.91 eV, L02 = 4.40 eV, T c = 308.3 K, and /3 = 0.142. 



III. MATERIAL PROPERTIES 

Numerical evaluation of the vdW energy and forces 
requires dielectric functions of the interacting materi- 
als. Accurate knowledge of the entire tensor i(cu, T) is of 
paramount importance. Before solving the electrostatic 
surface mode problem, material properties are discussed 
in this Section. 

For a growing number of substances, the dielectric 
functions on the imaginary axis are obtained through 
a Kramers-Kronig transformation of the absorption or 
reflection data followed by a fit to a multiple Lorentz 
oscillator models In this paper, the parameters given 
by Parsegiani 3 - and van Zwol and Palasantzas^i are 
used to construct dielectric functions of the following 
isotropic materials: silica^ (set 1); polytetrafluoroethy- 
lene^i (Teflon or PTFE); polystyrene^ (set 2), micai 3 - 
(Table L2.7, set b), gold 13 (Table L2.4, set 3), silveri£ 
(Table L2.5, set 1), and copperi 3 - (Table L2.6). 

There is much less information on e(i£, T) of liquid 
crystals. The usual difficulty of knowing the optical spec- 
tra in a wide energy interval is multiplied here by the 
need to know them separately for two polarizations and 
at different temperatures. Only a few liquid crystals have 
been studied experimentally well enough to enable a full 
model. One of the most studied compounds is 4-cyano- 
4-n-pentylbiphenyl (5CB), which will be used here as an 
exemplary positive uniaxial material. 

The dielectric functions of 5CB used in this paper 
are based on the three-band dispersion model developed 
by Wu and co-workers 1 20 > 22 n S4 The model accurately de- 
scribes the experimentally measured refractive indices in 
the (0.4-0.8) /im spectral interval for the entire temper- 
ature interval of the nematic phase 295.3 K < T < 308.3 
K. According to the model, the index dispersion in the 
visible is governed by three electronic transitions: one 
(7 — y (j* transition with Ao ~ 0.135 /im, and two ir — > it* 
transitions with Ai = 0.210 /<m and A2 = 0.282 /im. 
The oscillator strength of the a — > a* transition is very 
weakly temperature dependent. It can be extracted from 
the dispersion of the isotropic part of the refractive in- 
dex. The oscillator strengths of the ir — >■ ir* transitions 
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FIG. 2. Dynamic dielectric functions of liquid crystal com- 
pound 5CB, Eqs. ©-(O for the three temperatures 298.2 K, 
305.7 K, and 307.9 K. e" decreases with temperature, while 
e x increases with temperature. 



are temperature-dependent. In particular, they define 
the temperature variation of birefringence. Converting to 
imaginary frequency and using e^ ,J -(i^,T) = n\ (i£,T), 
the model reads 
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Here w = 9.19 eV, wi = 5.91 eV, and oj 2 = 4.40 
eV. The refractive indices of 5CB were measured by 
polarized UV spectroscopy^ and tabulated by Wu et 
al.— The temperature-dependent coefficients C extracted 
from those data are given in Table fl] 

Rotational relaxation of 5CB molecules and other low 
frequency processes are neglected here based on the fa- 
miliar argument 19 that their characteristic energies of less 
than 0.01 eV are much smaller than the first bosonic fre- 
quency at room temperature £1 « 0.16 eV. The conclu- 
sion is Eqs. ([H])-© represent the entire dynamical part 
of the dielectric functions. The functions are plotted in 
Fig. 

Static response requires a separate treatment. The 
static dielectric constant can be split into a temperature- 
independent isotropic part (w 10.7 for 5CB) and a 
temperature-dependent birefringent part. It is reason- 
able to assume that the temperature dependence comes 
from the order parameter. According to Li and Wu^. 
the order parameter of 5CB follows a universal relation 
oc (1 — T/T C )P, where T c = 308.3 is the nematic-isotropic 
transition temperature and (3 — 0.142 is a universal expo- 
nent. Adjusting the dielectric constants to the measured 
experimental values are lower temperatures^ one arrives 
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FIG. 3. Static dielectric functions of liquid crystal 5CB, 
Eqs. ©-©■ 



at the following parameterization 



on either side of one interface, and the full amplitude is 
found by matrix multiplications. This method is partic- 
ularly suited to the geometry at hand with its four in- 
terfaces. In the rest of this Section, the two fundamental 
transfer matrices are derived. 

The general solution in each region is a linear combi- 
nation of the transmitted and reflected waves. Consider 
interface (21). The solution in medium region 1 is sought 
in the form 



( 



Aie e i + Si e e i 
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J(q x x+q v y) 
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Here A± and B\ are unknown amplitudes, ef z in the ex- 
ponent is introduced for convenience, and the last ex- 
ponential factor reflects the uniformity of the Laplace 
equation in the (xy) plane. The exponents follow 
from Eqs. © and ©: 



q y ef z )±pi , 



(10) 



4cb(0>T) = 10.7 + 13.0(1 
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4cb(0,T) = 10.7- 6.5(1 



T 
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with T c = 308.3 K and f3 — 0.142. These functions are 
plotted in Fig. [3j 

To close this Section, several comments are in order, 
(i) The dynamical model (JS])-© is defined only at 5 dis- 
crete temperatures listed in Table |TJ It seems reasonable 
to extend the model to any temperature by fitting the co- 
efficients to the same universal factor a + b(l — T '/T C )P 22. 
This is not done in the present work, (ii) It is possible to 
convert Eqs. ([5|)-([6]) to a more familiar form by expanding 
the square and refitting the function to a linear combina- 
tion of oscillators. Such an additional fitting procedure 
might introduce undesirable errors, and therefore it is not 
employed here, (iii) Wu et al2£ provided a data set for 
another liquid crystal compound 5PCH, thus enabling a 
similar three-band dielectric model. 



IV. TRANSFER MATRICES 

In this Section, the anisotropic Laplace equation @ is 
solved. In analogy with the optics of multilayered media, 
it is convenient to construct solutions out of individual 
transfer matrices. Each matrix links wave amplitudes 



Pi 



1x ( e i e i) (cos 2 Q\ sin 2 "01 + co s 2 V'l) 
+ 9 2 (ei) 2 (sin 2 6i sin 2 ^i) 
+ 2q x q y £f (e\ — (sin 2 9\ cos-0i sin^i) 

+ q 2 y (er4) (cos 2 6»i cos 2 ipi + sin 2 ipi) 

+ g 2 (e r ) 2 (sin 2 1 cos 2 1 ) . (11) 

The quantity p is defined as the positive square root of p 2 . 
p is a function of the interface momentum components q x 
and q y , optical axis angles 9 and ip, imaginary frequency 
£, and temperature T. 

In slab region 2 the material is isotropic and the 
Laplace equation simply yields 



02 = (A 2 e qz + B 2 



-qz\ pi(<lxx+q y y) 



Q 



(12) 



(13) 



The matching conditions at z = £21 include the equal- 
ity of the transverse components of the electric field 
E x = —d(j)/dx and E v — —d(f>/dy (which lead to identi- 
cal relationships), and the equality of the normal compo- 
nents of D z = —e zk (dcf)/dxk)- The resulting two equa- 
tions can be rearranged to express the wave amplitudes 
of region 2 via the wave amplitudes of region 1: 
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The last equality defines the transfer matrix M 2 \ between the medium region 1 and slab region 2. A similar transfer 
matrix defines scattering at the (34) interface after index substitution 1 — > 3 and 2 — > 4. On the other hand, at the 
interface (23), the waves are incident from slab region 2. From the matching conditions one expresses the medium 
amplitudes A 3 and B 3 via the slab amplitudes A 2 and B 2 . After some algebra one obtains 
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A similar matrix describes scattering at the (54) inter- 
face. 

V. ONE SLAB 

In an isotropic liquid, a parallel-plate slab does not 
experience any macroscopic forces or torques. In an 
anisotropic liquid, the dependence of the dispersion forces 
on the inclination angle will result in a torque acting on 
the slab. For planar and other non-homeotropic surface 
alignments, the orientation of the optical axes on the two 
sides of the slab may in principle be different and this ef- 
fect also warrants analysis. If, in addition, the solid mate- 
rial itself is anisotropic, there will be another torque that 
will rotate the plate around its normal. The latter effect 
is not considered in this paper. (Note that all the cases 
mentioned are different from the mutual torque between 
two anisotropic bodies studied by Kats^ and Parsegian 
and Weiss J^) 

To determine the vdW energy of a single slab in an 
anisotropic host, only three regions of Fig. [T] need to be 
taken into account, for instance 1, 2, and 3. The scat- 
tering problem involves two transfer matrices 

(£)=#»-tf»(£)- (16) 

Surface states are defined as exponentially decaying at 
infinity. Accordingly, the amplitudes A± and B 3 must be 
set to zero. The top equation of Eq. (jTHJ) links the wave 
amplitudes A3 and B\ on either side of the system and 
hence defines the spatial structure of the surface mode. 
The bottom equation has the form W ■ B± =0. For a 
non- vanishing Bx, this implies W — 0, which yields the 
surface mode spectrum. Expressed via matrix elements 
of the transfer matrices, the spectrum equation is 

MllMH + MllMH = . (17) 



Substituting the matrix elements from Eqs. (fT4")) and (fT5)) 

it becomes 

W* = l- (g e '-**)(g e =»-Pi) . e -* = . (18) 
(qe 2 +P3){qe 2 + Pi) 

The free energy is then obtained as follows: (i) The spec- 
trum equation (fl8|) is substituted in Eq. ([3]); (ii) Polar 
coordinates q x — qcosx, q y = <?sinx, are employed in 
the integral over the surface vector; (iii) A new function 
u r = Pr/q is introduced. It is a function of the momen- 
tum angle \ but not of the momentum amplitude q. The 
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FIG. 4. Tilt Hamaker constant (|21[) of a parallel-plate solid 
slab immersed in 5CB, as a function of the optical axis tilt an- 
gle 9. The tilt is the same on both sides, 9i =63. The azimuth 
angles are ipi = ip3 = 0. The Rs are referenced from their 
respective values R(0) = -14.40, -12.90, -9.21, -4.25 zj 
for silica, Teflon, polystyrene, and mica, respectively. Teflon 
favors homeotropic alignment, 9 = 0, while other materials 
favor the planar alignment, 9 = n/2. Absolute temperature 
is T = 298.2 K. 
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FIG. 5. Same as Fig.[4]but for several metals in 5CB. The ref- 
erence tilt Hamaker constants are R(0) — —103.63, —116.13, 
—244.62 zj for gold, silver, and copper, respectively. All met- 
als favor the planar alignment, 6 = n/2. 



explicit form of u r follows from Eq. (jlip . (iv) The log- 
arithm is expanded in an infinite scries and integration 
over q is performed. The resulting expression for the free 
energy per unit area is 



F{- 
where 



kT V 2 - 1 f 27I dx / a . . 

2 n= Q m =l Jl > 
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(19) 
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and the same formula holds for A23 with 3 replaced by 
1. By analogy with vdW interaction between two semi- 
infinite bodies, the overall 1/d 2 dependence can be iso- 
lated by introducing a Hamaker-like constant R defined 
as F[ = i?/(127rd|). One arrives at 

n=0 m=l J ° 

(21) 

Note that R is defined with a negative overall sign to re- 
tain the visual appeal of energy profiles: large negative R 
implies lower energy and preferred orientation. In the fol- 
lowing, R will be referred to as "tilt Hamaker constant" 
to reflect its relation to the inclination of the optical axis. 

Figure 0] shows tilt Hamaker constants for several di- 
electric materials immersed in 5CB. The tilt angle is the 
same on both sides of the slab, Q\ = #3 and ipi = ipa = 0. 
For better presentation, the constants are referenced 
from their respective values at homeotropic alignment 
9 = 0. The reference values are listed in the cap- 
tion. Among the materials studied, Teflon favors the 
homeotropic alignment, while all other materials favor 
the planar alignment 6 — n/2. Similar plots for gold, 



FIG. 6. Tilt Hamaker constant (|21|) for single slabs in 5CB. 
The alignment is planar on both sides of the slab, 9i = 9z = 
tt/2 but the azimuth angle difference is systematically var- 
ied. The reference values are R(n/2,0) = -14.93, -11.92, 
-9.90, -5.06, -108.18, -119.87, -250.58 zj for silica, Teflon, 
polystyrene, mica, gold, silver, and copper, respectively. All 
materials favor the parallel alignment of optical axes, ip\ = z/>3 
and do not favor the cross alignment 163— ipi = tt/2. Absolute 
temperature is T = 298.2 K. 



silver, and copper are presented in Fig. [5j All the met- 
als favor the planar alignment. The absolute difference 
between the planar and homeotropic orientations is 4-6 
zj, i.e. almost an order of magnitude larger than for the 
dielectrics. 

For the planar surface alignment (as well as for any 
nonzero tilt angle), the azimuth orientation of the op- 
tical axes on the opposite sides of the slab can be dif- 
ferent. It is of interest therefore to investigate the vdW 
free energy as a function of the azimuth misalignment 
ip = ip3 — ipi- Such dependencies are shown in Fig. [B] 
for the planar alignment B\ = #3 = tt/2. For all di- 
electrics and metals, the parallel orientation of the op- 
tical axes, tpi = t/>3, is energetically preferred to the 
cross-orientation ^3 — tp\ — tt/2. However, this effect 
is relatively small, on the order of 0.1 zj. Azimuth mis- 
alignment makes a small contribution of to the overall 
energy balance. It is probably smaller than the variation 
from the uncertainty in the material parameters. 



VI. TWO PARALLEL SLABS 

In this Section, two parallel slabs immersed in a liquid 
crystal are considered. First, the spectrum of surfaces 
modes is derived from a product of four transfer matrices. 
The cases of semi-infinite slabs and finite-thickness slabs 
are analyzed in order. 
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A. 



General expression for the van der Waals energy expression for A is 



The system consists of five spatial regions separated 
by four interfaces, cf. Fig. [T] Collecting scattering at 
all four interfaces, the wave amplitudes in region 5 are 
expressed via the wave amplitudes in region 1 as follows 
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Surface states are defined by setting A\ = B 5 = 0. The 
top equation of Eq. (|2"2"|) defines the spatial structure of 
the surface mode. The bottom equation defines the spec- 
trum. Developing the bottom equation via matrix ele- 
ments one obtains 

{Ml\Mll + M 22 M 21 3 ) (MllNlH + M£mH) + 
(Ml\Mll + M 22 M 22 ) {MllMll + M 22 M 22 ) = 0. (23) 

Substituting here the explicit matrix elements from 
Eq. (|T4|) and (fT5j) and cancelling common positive- 
definite factors [this does not affect the final force after 
taking the logarithm in Eq. (O], the spectrum equation 
becomes 

W = [l-A 23 A 2ie - 2 ^][l-A43A45e- 2 ^] 

- cT 2 ^ [A 23 - A 21 e- 2qd *] [A 43 - A 45 e- 2 ^].(24) 

Here L = L32 — L43 is the gap between the slabs and the 
factors A 43 and A45 are defined according to Eq. (|2"0)l 
with 2 replaced by 4. 

In accordance with the recipe ([3]), the free energy per 
unit interface area is 



This formula can be used for numerical calculations. The 
input parameters are the geometrical factors (I2/L, di/L, 
the orientation of optical axes 9 r and tp r and the temper- 
ature T. The last three parameters define the quantities 
u T that enter via A^-. 

In the limit of thin slabs the integrand is nonzero 
within the large interval < t < L/d so the second term 
under the logarithm does not contribute much. Then A 
converges to the quantity R defined in Section [V] mul- 
tiplied by the factor (L/d) 2 that is responsible for the 
difference in definitions of A and R. 



B. Two semi-infinite slabs 

Two semi-infinite bodies interacting via a gap L is 
the basic vdW geometry. In this Section, the Hamaker 
constant for uniaxial anisotropic media is derived from 
the general formalism and then numerical results are 
presented. The spectrum of surface modes is given by 
Eq. (|2"o]) . Going over to polar coordinates, expanding 



(2tt) 2 L J 



(25) 



If the gap between the slabs is large, L 3> ^2,^4, then 
the second term in Eq. ([24]) vanishes. The first term 
under the logarithm in Eq. (|25[) splits into two parts, 
each corresponding to the vdW energy of an isolated slab 
surrounded by the medium. The total energy reduces to 
a sum of two terms derived in the preceding Section. 

If the slabs are thick, di^d^ 3> £, the spectrum equa- 
tion (l24l) reduces to 



W c 



1 - A 23 A 43 e c s~ 



. 



(26) 



It will be analyzed in Section IVIBI Here the general ex- 
pression (|25p is adapted for numerical evaluation. 

In the integral over q, polar coordinates q x — qcosx, 
q y = q sin \ are useful. Then the Hamaker "constant" 
A = 12itL 2 Fi is introduced to account for the usual 1/L 2 
scaling of the energy. It is also convenient to change 
the integration variable from q to t = qL. The final 
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FIG. 7. Hamaker constant (|28[) for two semi-infinite bod- 
ies made of different dielectric materials separated by liq- 
uid crystal 5CB at T = 298.2 K. The reference values are 
A°°(0) = -16.09, -16.57, -11.87, and -5.32 zj for silica, 
Teflon, polystyrene, and mica, respectively. The attraction is 
strongest at 63 — 0. 
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FIG. 8. Hamaker constant (|28(l for two semi-infinite bodies 
made of gold, silver, and copper separated by 5CB at T = 
298.2 K. The reference values are A°°(0) = -128.82, -144.94, 
and —283.55 zj, respectively. The attraction is strongest at 
6*3 = for all metals. 



the logarithm and integrating over q results in 

r2 "d X (A 23 A 43 ) r 



3kT 



oo 
2 \ " 



1 



2n 



Li lit lr\ / L !.!• > 

(28) 

Figure [7] shows the Hamaker constant A 00 for several 
dielectric materials separated by 5CB, as a function of the 
tilt angle of liquid crystal's optical axis. For the conve- 
nience of presentation, the Hamaker constants have been 
referenced from their values at 9$ = 0. The reference 
values are listed in the caption. All the materials show 
preference of homeotropic alignment #3 = 0. Among the 
materials studied, Teflon has shown the largest differ- 
ence between the homeotropic and planar vdW energies 
(about 5.5 zJ). 

Metals possess qualitatively similar angle variations 
of vdW energy, as shown in Fig. [8l The attraction is 
stronger for the homeotropic alignment. However, the 
absolute scale of the variation is about one order of mag- 
nitude larger than for dielectric materials. A stronger 
vdW attraction for homeotropic alignment may be a gen- 
eral feature of 5CB and perhaps of any positive liquid 
crystal. 

One expects the inclination dependence to go away as 
the temperature increases and the medium becomes opti- 
cally isotropic. The temperature dependence of A 00 (#3) 
for the silica-5CB system is shown in Fig. [§] 



C. Two finite-thickness slabs 
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FIG. 9. Hamaker constant (|28[) f° r two semi-infinite silica 
bodies separated by liquid crystal 5CB for five temperatures. 



vdW free energy. At large separations L ^> d 2 ,d,4, one 
expects one-slab effects to be dominant. For all ma- 
terials except Teflon the vdW energy is minimal when 
the liquid crystal is planar aligned on both sides of the 
slab. In the opposite limit of very small separations, 
L -C d2,d 4 , one expects the interaction effects to dom- 
inate. As discussed above the interaction energy favors 
the homeotropic alignment for all the materials studied. 
It is of interest therefore to study the evolution of the 
optimal orientation of the middle liquid crystal and to 
follow the interplay between the single-slab and interac- 
tion effects. 

Consider the 5CB-silica-5CB-silica-5CB system as an 
example. A single silica slab in 5CB favors the planar 
alignment on both surfaces. Accordingly, one sets #1 = 
f? 5 = 7r/2, and 4>i = ip3 = ^3 = 0. The tilt angle of 
the middle LC section 6*3 remains variable to include the 
possibility of homeotropic and other alignments. The 
slabs are assumed to be of the same thickness, d 2 = d 4 = 
d. The overall vdW energy is studied as a function of 6*3 
for different ratios d/L. 

Results of numerical calculations are shown in Fig. [TOl 
For large d/L > 5, the vdW energy is dominated by 
the interaction across the gap, and the f9 3 dependence is 
virtually the same as in the semi-infinite case. (Compare 
the d/L = 50 plot in Fig. [TO] with the T = 298.2 K plot 
in Fig. [HI) As the slabs get thinner, the single-slab effects 
grow stronger and eventually dominate. Accordingly, the 
planar orientation 6*3 — tt/2 becomes the absolute energy 
minimum at d/L < 0.5. At intermediate thicknesses, 
0.5 < d/L < 1.0, the energy has two local minima, at 
6*3 = and 83 = tt/2, as can be seen in the figure. 



To analyze a system of two slabs oriented parallel to 
each other, the full four-transfer matrix solution (|27l) is 
needed. The present study is focused on finding an op- 
timal orientation of the optical axes that minimizes the 



VII. SUMMARY AND DISCUSSION 

Collective behavior of colloidal particles in anisotropic 
media is a fascinating and complex subject. This be- 
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FIG. 10. Hamaker "constant" A, Eq. (|27[) . of two silica slabs 
in 5CB for several slab thicknesses. The alignment of liquid 
crystal on the outside surfaces is planar, 9i = 65 — it/2, while 
the alignment in the gap is varied. The reference values are 
-16.09, -39.21, -54.25, -67.33, -87.77, -122.12, -186.06 
zj, (for d/L going from large to small). The absolute temper- 
ature is T = 298.2 K. 



havior is determined by the balance of surface alignment 
energy, bulk elastic energy, electrostatic forces, van der 
Waals forces, and thermal fluctuations. Given the tech- 
nological importance of both liquid crystals and non- 
aqueous colloids it would be desirable to have a com- 
prehensive theory of colloidal stability in liquid crystals 
of the same clarity as the classical theory of colloidal 
stability in isotropic fluids^ Such a theory does not yet 
exist. 

The main purpose of the present work has been to 
demonstrate that even a single component in this mix, 
the van der Waals force, is complex and can lead to non- 
trivial consequences. Due to the complexity of the gen- 
eral problem, only the plane geometry in the nonretarded 
limit has been analyzed. The vdW free energy can be ob- 
tained in this case from the spectrum of electromagnetic 
surface modes relatively easily. Unlike previous works, 
the focus here has been the vdW energy dependence on 
the inclination angle of the optical axis. 

A significant barrier for any realistic calculation of 
vdW forces is the lack of reliable parameterization of the 
dielectric function on the entire imaginary frequency axis. 
In liquid crystals, this is further complicated by birefrin- 
gence and a strong temperature dependence. In this pa- 
per, a three-oscillator temperature-dependent model of 
5CB has been constructed based on the real-frequency 
data of Wu and co-workersi 20 i 22 ~ — More work will be 
needed to further validate and refine the model presented 
in Section Hill 

In an anisotropic fluid, the vdW energy of a parallel- 
plate slab becomes a function of the tilt angle between 
the optical axis and the surface normal. Energy profiles 
have been calculated in Section [V] All studied materials, 



except Teflon, favor the planar alignment of the optical 
axis. If the real anchoring orientation is different from 
the optimal one, the dependence of the vdW energy on 
the tilt angle will generate a vdW torque that needs to be 
taken into consideration in determining the equilibrium 
orientation of the slab. The torque disappears as the 
temperature is raised above the nematic-isotropic transi- 
tion. 

In the case of planar alignment and other non-homeo- 
tropic surface alignments, the optical axes may have dif- 
ferent azimuth orientation on either side of the slab. The 
vdW energy is in general a function of the azimuth angle 
mismatch. This effect has been investigated in Section fVl 
and found to be numerically small. All the materials 
studied favor parallel orientation of the optical axes, i.e. 
equal azimuth angles on both surfaces. 

When two slabs are brought close together, they at- 
tract via a vdW force that is a function of the optical 
axis direction of the intermediate medium. A general so- 
lution to this problem has been developed in this paper, 
cf. Section l"VI A[ It has been found that the vdW force 
is strongest for the homeotropic orientation of 5CB for 
all the materials analyzed. 5CB is a positive liquid crys- 
tal. Thus the main result suggests the vdW attraction is 
strongest when the surface normal is parallel to the line 
of largest polarizability of the medium molecules. 

The last observation might have important implica- 
tions for the stability of colloids in liquid crystals. Imag- 
ine a pair of spherical particles in a positive liquid crystal 
under weak anchoring conditions. When the center-to- 
center line is parallel to the optical axis the vdW attrac- 
tion will be stronger than when it is perpendicular to 
the optical axis. Since the electrostatic repulsion from 
double-layer overlap also depends on the dielectric con- 
stant, the balance between the attractive and repulsive 
forces will depend on the mutual orientation of the par- 
ticles. As a result, the particles may attract along some 
directions but repel along others, which could lead to 
chain formation. This intriguing possibility warrants fur- 
ther investigation. 

Finally, the effects of finite slab thickness have been in- 
vestigated. Using the general solution, Eq. (P2T1) . smooth 
evolution of the vdW energy from the gap dominated 
limit to the slab-thickness dominated limit has been ob- 
served. At least for some materials (such as silica in 
5CB) this implies that the planar alignment in the gap 
is preferred for thin slabs and large gaps, while the 
homeotropic alignment is preferred for thick slabs and 
small gaps. 
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